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Abstract 

Using the Laplace sum rule (LSR) approach, which is less affected by the contribution of 
the higher mass hadronic states than the Finite Energy Sum Rule (FESR), we test the 
reliability of the existing estimate of the K°-7C° mixing parameter from the four-quark 
two-point correlator. We obtain, for the renormalization group invariant 5-parameter 

-i 2 - 

fK/(l-2fn) Bk, the upper bound: 0.83 and the best estimate: 0.55 ± 0.09. Combining 
the previous estimate with the updated value of JbVBb = (1-49 ± 0.14)/,,- obtained from 
the same LSR method, one can deduce the best fitted values (p, rj) « (0.09,0.41) of the 
CKM parameters. 
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1 Introduction 



A reliable estimate of the Bx parameter is of prime importance for understanding CP- 
violation. In a series of papers Pich and de Rafael and coworkers [[!], 0] have estimated its 
value using the Finite Energy Sum Rule (FESR), from the two-point correlator: 

M>(g 2 ) = • J d 4 xe^(0\TO AS =2(x)Oi s=2 (0)\0), (1) 
associated to the AS = 2 four-quark operator: 

O as =2 = s^Ldsj^Ls, (2) 
where L = (1 — 7s)/2, while the -B^-parameter is defined as: 

(K°\O A s = 2\K ) = p 2 K M K B K , (3) 

where fx — 1-2/^ is the kaon decay constant (f n = 93.3 MeV). The renormalization group 
invariant (RGI) 5-parameter is defined as [l|]: 

B K = B K (a s (/i 2 ))~ 2/9 , (4) 

and its latest improved value from FESR is 0: 

B K = 0.39 ±0.10 (5) 

for fx — 1-28/,,- MeV. The improvement with respect to the original paper in resides 
mainly in the hadronic parametrization of the spectral function where, within the frame- 
work of chiral Lagrangian and Weinberg sum rules, higher mass-state contributions have 
been included into the spectral function. Indeed, in the FESR approach, their role is im- 
portant as this sum rule is sensitive to the high-energy behaviour of the spectral function. 
A priori, this is not the case of the Laplace sum rule (LSR), as the exponential factor 
entering in it enhances (when it is operative) the low-energy contribution of the spectral 
function. For this reason, we might expect that the Laplace sum rule will be less sensitive 
to the higher mass hadronic state effects, which are not under good control. 



2 The four-quark two-point correlator 

The QCD expression of the spectral function has been evaluated in |l| including dimension- 
four operators, while the radiative perturbative correction has been re-evaluated recently 
0. The RGI two-point correlator reads: 

1 "nnn, s 1 14,, , f 3649 a s m 2 Q Vt 1 

_ Im *«c D(() = {aM) V, | t _ _^ _ 40 _, _ _ + ... I ? (6) 

where: 

Q, = 200m* - 207r[167rm s (ss) - (a s G 2 )]. (7) 
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We shall use the following values of the QCD parameters: 

A 3 = (0.35 ±0.10) GeV, (8) 

corresponding to a s (M z ) = 0.118 ± 0.006 to two-loop accuracy. We shall use j|: 

/ 2tt \ ~ 4/9 / a \ 

m s = (270 ± 18) MeV m s = m s [ 1 + 0.89— + ... (9) 

V9a s / V 7r / 

K ' 0.6-1 (dd) = — (188.3 MeV) 3 ( ] , (10) 



(dd) ' ' \9a 

and: 

(a s G 2 ) = (0.06 ± 0.03) GeV 4 . (11 

The hadronic part of the spectral function will be parametrized in the same way as in 
which is given in Eq. (25) of that paper (after correcting obvious misprints): 

ii W 11 288tt 2 \f n J Jt w Jt 20 V ' t ' t 

t-i + h-A 2 Iimn^rto-imn^fe) 

t t ) IX 71 

+ 2\(l, ^, ^) -ImU^^-lmU^ih) 

V t t J 7T 71 



+ 



tl t2 A 2 



J + J~ 1 ) + t 



-imnw^o-imn^^) 

71 7T 



+ &(t-t c )-lm^ QCD (t), (12) 

7T 

where the index 2=0, 1 refers to the hadronic states with spin 0, 1, and: 

-Imn^(t) = -lmU { y\t) + -ImU { 2(t), (13) 

7T 7T 71 

where the index V, A refers to the vector and axial-vector currents; A 1 / 2 is the usual 
phase space factor. We shall use the parametrization with finite width corrections (FWC) 
proposed in Ref. 0], but we shall compare the results with the ones obtained in the 
narrow width approximation (NWA). 



3 The Laplace sum rule analysis of Bk 

We shall be concerned with the Laplace sum rule: 

£(r) = dte- tT -lm^/(t), (14) 
Jto vr 

where to is the threshold specific for each channel, while t c is the one of the QCD contin- 
uum coming from the QCD spectral function given previously. The sum rule analysis is 
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standard. We are looking for the region of stability^ (minimum sensitivity)^, in the sum 
rule variable r and study the independence of this r-stability (stability in t c corresponds 
to the minimum sensitivity in the higher meson mass contributions). For illustrating the 
analysis, we show in Fig. 1 the r-stability for different choices of the QCD parameters and 
for various values of t c . It can be noticed that the stability inflexion point is at r = 0.125 
GeV -2 , which is exceptionally small compared with the hadronic scale 1 GeV -2 . This 
feature indicates that the parametrization of the spectral function by the KK alone will 
be a bad approximation, as expected from the FESR analysis. This is essentially due to 
the quartic t-behaviour of the QCD spectral function and to the huge non-perturbative 
corrections, which push to work at smaller r-values for the Operator Product Expansion 
(OPE) to make sense, but at this r-values, one is quite far away from the KK threshold. 
However, at this small r-value, a possible large, but quantitatively less controllable, as 
instanton-like effect appearing in high powers of r is obviously negligible. The result is 
shown in Fig. 2a (curve 0) where the value of the .B-parameter increases brutally with t c . 
The situation is already improved when one includes the KK* (curve 1), where a stability 
in t c occurs at the inflexion point. We have used the parametrization of the form factor 
including non-resonant Kr] states proposed in g], which we have compared with the one 
within NWA using the /^-coupling 7#* = (2.63 ±0.16) from r-decay data and from QCD 
spectral sum rule (QSSR)||, with the normalization: 

hmTl K *(t) = ^6(t-M^). (15) 

The FWC is almost negligible although it tends to pronounce the existence of the inflexion 
point. Then, the positivity of the spectral function implies an upper bound: 

f K / {1.2 f n )] 2 B < 0.83. (16) 

We add to the former the K*K* and KK\ effects (curve 2). We use the i^i-mass of 
1.39 GeV, while we determine its coupling 7^ ~ 4.6 using the first Weinberg sum rule 
involving the K\ K, K x and K* treated in a NWA. We use f K * ~ (28.3 ±4.5) MeV (recall 
/tt = 93.3 MeV) from QCD spectral sum rules (QSSR)[|J. For the curve 3 (analogous to 
curve 3 of 0), we add the effects of KKq and K*K\. Although there is a big change 
from curve 1 to 2 (a similar feature is obtained from FESR in Fig. 2 of @), there is 
a negligible effect from curve 2 to 3 (this is not the case of FESR). Here, we see that 
the exponential factor starts to be operational as the effects of the higher states become 
suppressed. Indeed, the inclusion of the effects of the K*Kq, K*K±, KqKq and K\K\, 
and some other combinations involving higher mass states give negligible effect (at the 
optimization scale, the exponential introduces a weight factor of about 0.3 instead of 1 for 
the FESR). Therefore, we include them into the QCD continuum. The previous feature is 
one big advantage of the Laplace sum rule compared with the FESR. Also by comparing 
curve 2 here and in Fig. 2 of [fj], we can notice that the result from the LSR is more 
stable in t c than the one from FESR. 

lr The meaning of stability has been discussed in Q by comparing the exact and approximate solutions 
of the sum rule for the harmonic oscillator in quantum mechanics. What happens in QCD is very similar 
to this case. 

2 At the minima or at the inflexion point, we have a balance between the continuum and the nonper- 
turbative contributions. 
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In Figs. 2, we have an inflexion point at the value of t c = 7.5 — 8 GeV 2 , which is inside 
the duality region 6.5-9.5 GeV 2 expected from the FESR JI], |2|]. We show in Fig. 2b the 
effect of the subtraction point where the QCD parameters are evaluated as well as the 
one due to the S'[/(3)-breaking of the quark condensate. Taking into account different 
sources of uncertainties, we deduce at the inflexion point: 



f K /{l.2U) B K = 0.55 ± 0.02 ± 0.05 ± 0.05 ± 0.05, 



(17) 



where the errors are due respectively to the choice of the subtraction point, the value of 
a s or A, the SU(3)-breaking of the condensate, and the localization of the inflexion point. 
Adding these errors quadratically, we finally obtain the best estimate: 



'B 



K 



0.55 ±0.09. 



A conservative estimate could be obtained inside the whole range of the t c -duality region 
from 6.5 to 9.5 GeV 2 . In this case, one obtains from Fig. 2b: 



fi</(l.2f n ) 



'B 



K 



0.58 ±0.22, 



(19) 



which, however, we expect to be a number with a too generous error. On the other hand, 
by combining the results in Eqs. (18) and (19) with the one in Eq. (5) (after rescaling 
the value of fx) from the FESR, we deduce the (weighted) average: 



f K /(1.2U) 



~B 



K 



0.51 ±0.07. 



(20) 



We have not included in this average the existing estimates from 3-point function sum rules 
(see e.g. ||] for a critical review on these estimates). Indeed, due to the more complicated 
analytic property (in particular, we loose the positivity property of the spectral function) 
of this quantity, the inclusion of the higher resonance mass states into the corresponding 
spectral function is therefore highly non-trivial (interference terms) and can even become 
uncontrolled. In this respect, we expect that the estimate from the 3-point function is 
less accurate than from the 2-point function. 



4 Conclusions and values of the CKM parameters 
from the LSR 

We have derived the value of the invariant I?-parameter using the LSR. We have obtained 
the upper bound in Eq. (16), the best estimate in Eq. (18) and the conservative estimate 
in Eq. (19). Our result is slightly higher than the one from FESR, but agrees with it within 
the errors. By combining the LSR and FESR estimates, we deduce in Eq. (20) the average 
from QSSR. Our previous result agrees with the estimate from the chiral Lagrangian ||. 
Our bound excludes some other estimates in the existing literature and raises again the 
question on the actual reliability of the value 0.8 from present lattice calculations J7J, 
which is, however, obtained in the SU(3)f symmetry limit. In our approach, the SU(3)f 
symmetry breaking is mainly reflected in the deviation of fx from f n (see Eq. (12)). 
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Then, if we ignore the SU(3)f breaking, i.e taking f K = f n , we would obtain B K ~ 0.79, 
in (surprising) good agreement with the present lattice value. 

Our estimate of Bk combined with the recent LSR estimates of fs |§, -Bb@ and the value 
of the top mass will certainly give stronger constraints on the CP-violation parameters 
That can be explicitly seen from various fits || [TOJ] . Indeed, by using our best estimate 
in Eq. (18) and our best value 0: 

f B \[B~B = (1.49 ±0.14)/,, (21) 

which are, actually, both obtained from only one method (LSR for the two-point func- 
tions), and using the standard Wolfenstein parametrization, one obtains for the CKM 
parameters of the unitarity triangle, the best fit for x 2 — 0-58 [O: 



p « 0.09 r, « 0.41, (22) 

and, in Fig. 3, the present allowed domain of the set (p, n). We can consider the previous 
set (p, rj) values as the first consistent (in the sense of a single type of sum rules used) 
estimate from QSSR. One can also notice from [ID| that p is more sensitive to the changes 
of the QCD input parameters than rj. 



Figure captions 

Fig. 1: r-stability of the invariant parameter Bk for different choices of the QCD 
parameters: a) running mass and coupling evaluated at £; b) the same as (a) but at 
r-variable; c) the same as (a) but at r = 0.125 GeV~ 2 ; d) the same as (b) but using 
(ss)/(dd) = 0.6. 

Fig. 2a: testability for different parametrizations of the spectral function: curve 0: 
KK; curve 1 : (0) + KK*; curve 2 : (1) + KK X + K*K*; curve 3: (2) + K X K* + KK*. 



Fig. 2b: Effects of different QCD parameters on the t c -stability. The label of each curve 
is the same as the one in Fig. 1. 

Fig. 3: Allowed domain in the p-?7-plane from a simultaneous fit of the present data 
compiled in Table 1 of [T(J and of the QCD parameters given in Eqs. (18) and (21). The 



triangle represents the best fit. 
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